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ABSTRACT 

We present a spectro-astrometric survey of molecular gas in the inner regions of 16 protoplanetary 
disks using CRIRES, the high resolution infrared imaging spectrometer on the Very Large Telescope. 
Spectro-astrometry with CRIRES measures the spatial extent of line emission to sub-milliarcsecond 
precision, or < 0.2 AU at the distance of the observed targets. The sample consists of gas-rich disks 
surrounding stars with spectral types ranging from K to A. The properties of the spectro-astrometric 
signals divide the sources into two distinct phenomenological classes: one that shows clear Keplerian 
astrometric spectra, and one in which the astrometric signatures are dominated by gas with strong 
non-Keplerian (radial) motions. Similarly to the near-infrared continuum emission, as determined by 
interferometry, we find that the size of the CO line emitting region in the Keplerian sources obeys a 
size- luminosity relation as Rqo oc L^"^. The non-Keplerian spectro-astrometric signatures are likely 
indicative of the presence of wide-angle disk winds. The central feature of the winds is a strong 
sub-Keplerian velocity field due to conservation of angular momentum as the wind pressure drives 
the gas outwards. We construct a parametrized 2-dimensional disk+wind model that reproduces the 
observed characteristics the observed CO spectra and astrometry. The modeled winds indicate mass- 
loss rates of > 10~^^ — lO~^M0yr~^. We suggest a unifying model in which all disks have slow 
molecular winds, but where the magnitude of the mass- loss rate determines the degree to which the 
mid-infrared molecular lines are dominated by the wind relative to the Keplerian disk surface. 
Subject headings: techniques: image spectroscopy — protoplanetary disks — line: profiles — ISM: 
molecules 



1. INTRODUCTION 

Planets are believed to predominantly form in the inner 
regions of circumstellar disks surrounding young stars, 
the so-called planet-forming zone (PFZ). Consequently, 
searches for embryonic planetary systems are expected to 
be most fruitful if they are focused on the inner regions 
of planet-forming disks, whether using direct detection 
methods or indirect approaches such as imaging signa- 
tures of planet-disk interactions. While typical gas-rich 
disks around young stars may extend to 100 AU, or more 
in a few extreme cases, the outer reaches are too tenu- 
ous to form planetary cores within the lifetime of the 
disk, unless the disk is massive enough to become gravi- 
tationally unstable (Boss 1997; Marois et al. 2008). A few 
mature planets or substellar objects have been imaged at 
large radii around A stars, but the majority of the known 
planetary systems are believed to have formed within 20 
AU (Pollack et al. 1996), followed by radial migration in- 
wards due to planet-disk interactions during the gas-rich 
phase of the disk (Alibert et al. 2004), consistent with 
the currently known distribution of exo-planets detected 
with the radial velocity and transit methods. 



1.1. An observational challenge 

A primary difficulty in studying planet formation in 
progress is the small angular sizes of the PFZs. The 
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nearest protoplanetary disks are located, with a few no- 
table exceptions such as the TW Hya association and a 
few Herbig Ae stars, at distances in excess of 100 pc. 
Resolved images of typical PFZs must be obtained at a 
spatial resolution better than O'.'l. Further compounding 
the problem is that many of the main tracers of PFZs are 
found in the infrared wavelength range (2-200 /im) be- 
cause the relevant gas temperatures are in excess of 100 K 
- making ground-based observations challenging. Gen- 
erally, specialized instrumentation is needed to obtain 
the requisite spatial resolution. In recent years, ground- 
breaking progress has been made in infrared interferome- 
try of the innermost regions of protoplanetary disks from 
1-2 /im (e.g., Millan-Gabet et al. 1999; Eisner et al. 2005) 
and near 10 /im (van Boekel et al. 2005). However, inter- 
ferometric observations of gas have thus far been limited 
to single band or relatively low spectral resolving power 
and to young A and B stars. This has, for the most 
part, limited interferometry to dust, gas continuum and 
hydrogen recombination lines (Tatulli et al. 2007, 2008; 
Isella et al. 2008; Eisner et al. 2009). 

Our understanding of the structure and dynamics of 
inner protoplanetary disks and PFZs remains limited. 
While there appears to be a common end result of disk 
evolution - the formation of a planetary system, com- 
plete with gas giants and perhaps smaller rocky planets, 
as well as, presumably, planetesimals, comets and zodia- 
cal dust - the pathway is unclear. It is known that disks 
exist and that they carry most of the angular momen- 
tum of a young stellar system, as evidenced by resolved 
sub-millimeter imaging of their outer regions (Koerner 
et al. 1993; Mannings & Sargent 1997; Qi et al. 2008), 
and matching the momentum distribution of the solar 
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system and other planetary systems. Their spectral en- 
ergy distributions show us that the disks are gas-rich 
throughout; if they were not, the disks would be flat, not 
flared with essentially pressure-supported scale heights 
(Kenyon & Hartmann 1987). 

The presence of accretion shocks and jets tell us that 
the disks are actively accreting requiring the gas to be 
sufficiently viscious to allow angular momentum trans- 
port (Koenigl 1991). Limited disk lifetimes of about 
6 Myr have been inferred, albeit with significant varia- 
tion, demonstrating the need for efficient dispersal mech- 
anisms (Haisch et al. 2001). Disks do not extend all the 
way inward to the stellar surface, but exhibit a complex 
structure due to the subhmation of dust at high temper- 
atures leading to rapid opacity gradients (Natta et al. 
2001), and a coupling of the disk to the stellar magnetic 
field facilitates both accretion and mass loss (Koenigl 
1991; Shu et al. 1994). 

Some disks appear to have large excavated inner re- 
gions of low dust opacity (Strom et al. 1989), compared 
to the outer disk, but it appears that not all disks go 
through such a stage (Muzerolle et al. 2010), and it is 
still an open question whether this is directly related to 
the formation of planetesimals or even giant planets, or 
whether another dispersal mechanism is at play (Alexan- 
der & Armitage 2009). Clearly, observing what the gas 
actually does in the inner disk impacts our understand- 
ing of all these disk properties and enables us to esti- 
mate their relative importance for the evolution of the 
disk. Ultimately, the motion of the inner disk gas may re- 
veal the presence of accreting protoplanets (Regaly et al. 
2010). 

1.2. This paper 

We present a spectro-astrometric imaging survey of 
molecular gas in the PFZs of disks around a sample of 
solar type stars using CRIRES on the European Very 
Large Telescope (Kauff et al. 2004). The primary goal 
is to directly determine the basic distribution and kine- 
matics of the gas and to relate this to the process of 
planet formation and inner disk evolution. For instance, 
in a purely passive, non-accreting, disk it may be ex- 
pected that the gas orbits at essentially Keplerian speeds, 
dictated only by the mass of the central star, neglect- 
ing minor corrections for the mass of the disk itself 
{dV/VKepier ^ 10 ~^) and for the radial pressure gradi- 
ent in the disk {dV/VKepier ^ 10"^ Cuzzi et al. 1993). 
The presence of double-peaked emission line profiles from 
protoplanetary disks can indeed be explained by gas in 
Keplerian orbits (Carr et al. 1993; Blake & Boogert 2004; 
Pontoppidan et al. 2008). However, disks are in general 
not passive as they accrete and eventually dissipate, so 
departures from pure Keplerian motion are expected at 
some level, and is indicated in FU Ori stars for the CO 
bandhead at 2.3 /im (Calvet et al. 1991; Martin 1997; 
Hartmann et al. 2004). Do mid- infrared molecular lines 
- in fact - trace material strictly in Keplerian motion, or 
are there significant non-Keplerian components present, 
and, if so, are they dominated by infall or outflow mo- 
tions? What disk radii are traced by rovibrational CO 
lines; the inner edge of the disk at 0.01-0.1 AU, the ter- 
restrial planet region at 0.1-1.0 AU or the giant planet 
region at 1.0-10 AU? What is the origin of the absorp- 
tion components seen in CO rovibrational spectra - ab- 



sorption from edge-on disks, foreground clouds, remnant 
envelopes or disk winds? 

Our spectro-astrometric survey was primarily focused 
on what is probably the best tracer of molecular gas in 
the inner disk, or at least the most easily observable 
- the fundamental rovibrational band of CO centered 
at 4.67 /im. These CO fines are traditionally thought 
to trace warm gas in disks at ~1 AU, based on typical 
line widths and excitation temperatures (Najita et al. 
2003; Blake & Boogert 2004), and are invariably bright 
in nearly all classical T Tauri and Herbig Ae stars. One 
of the most basic outcomes of this survey is the direct 
measurements of the size of the CO line emitting regions. 

We demonstrate that astrometric signals in CO were 
detected for all sources on AU-scales, but with varying 
amplitude and with an intriguing range of structure: The 
CO spectra are divided into three rough phenomenolog- 
ical classes, based on the line shape in combination with 
the shape of astrometric spectra: Keplerian disks charac- 
terized by double-peaked line-proflles, single-peaked line 
sources with broad wings, and self-absorbed sources. 

This paper is arranged as follows: In §2 the survey and 
data reduction are described, including a detailed discus- 
sion of the capabilities of the spectro-astrometric mode 
of CRIRES for super-resolution imaging. In §3), we dis- 
cuss the results of fitting simple Keplerian disk models 
to the data. The central issue is that many CO line 
and astrometric spectra cannot be explained by Keple- 
rian velocity fields. We introduce a non-Keplerian class 
of emission lines in §4 and explain why a purely Keplerian 
model fails. In §5, we develop a 2-dimensional model that 
adds a disk wind to a regular Keplerian, fiared disk, and 
demonstrate that this provides a framework for match- 
ing all CO line and astrometric spectra from classical T 
Tauri stars under specific circumstances, namely if the 
wind is slow and uncollimated. We suggest that there is 
a smooth transition from lines dominated by Keplerian 
motions to wind-dominated lines, possibly scaling with 
the mass-loss/accretion rates. In §6 the implications for 
our understanding of disk dispersal and the nature of the 
warm molecular disk surface layer are discussed. 

2. OBSERVATIONS 

Spectro-astrometric observations were obtained as part 
of a large CRIRES survey of infrared molecular emis- 
sion from protoplanetary disks and young stellar objects 
within the framework of the European Southern Obser- 
vatory (ESO) Large Program 179.C-0151 (Pontoppidan 
et al. 2011). Spectro-astrometry is a highly sensitive 
method that allows a single telescope to obtain both spa- 
tial and kinematic information on gas-phase lines on very 
small scales, < 1 milliarcsecond, and at very high spec- 
tral resolution. A/ A A ~ 100 000. The final accuracy of 
a spectro-astrometric measurement depends linearly on 
both the signal-to-noise and the width of the spatial PSF. 
Basically, the method measures the spatial centroid off- 
set of the spectrum as a function of wavelength across a 
line or other spectral feature, relative to the continuum. 
This approach can reveal spatial structure on scales much 
smaller than the formal diffraction limit of the observa- 
tion. Spectro-astrometry was first developed for chro- 
matic imaging and spectroscopy for the detection of stel- 
lar binaries in the visible range using specialized instru- 
mentation (Beckers 1982; Christy et al. 1983; Aime et al. 
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1988). The modern form, using an echelle spectrograph, 
was first presented by Bailey (1998), and is reviewed by 
Whelan & Garcia (2008). Infrared (A > 1 /im) spectro- 
astrometry of molecular lines with CRIRES was intro- 
duced by Pontoppidan et al. (2008), who presented CO 
data from three transitional disks, TW Hya, HD 135344B 
and SR 21. They showed that sub-milliarcsecond preci- 
sions could routinely be achieved and that the basic ge- 
ometries of the line emitting regions - sizes, inclinations 
and position angles - could be determined with a high 
degree of confidence. Here we extend this sample to a 
much wider range of disks in terms of stellar type and 
evolutionary stage. 

2.1. Observing strategy 

Targets were selected for spectro-astrometric obser- 
vations according to overall brightness and line-to- 
continuum contrast, as well as to cover as wide a range 
as possible in known disk and stellar characteristics. The 
line- to- continuum contrast is a particularly important 
parameter to consider, since the accuracy of the mea- 
sured astrometric signal depends roughly linearly on this 
parameter (see §2.3). However, we were successful in ob- 
taining high-quality spectro-astrometry for sources with 
line-to-continuum ratios spanning {F — C)/C ^ 0.1 — 2, 
where F is the total flux and C is the continuum flux. 
The central stars include spectral types from K7 to late 
A, and cover luminosities of 0.1 to 100 L©. The tar- 
gets also span a wide range in accretion rates from 
logM ^ -9 to -6M0yr-^ 

Overall, the sample is intended to represent the di- 
versity found among protoplanetary disks to the extent 
that is possible given the limitations on the size of the 
sample. A log of the observations is given in Table 1, 
including those already described in Pontoppidan et al. 
(2008), while Table 2 summarizes the properties of the 
central stars. 

The lack of a correction for differential refraction be- 
tween the effective wavelengh of the slit viewing camera, 
which usually operates at H- or i^-band, and the spectral 
wavelength precluded the observation of low-elevation 
targets for much of duration of the campaign^, includ- 
ing sources in the Taurus and Chamaeleon star forming 
regions. As a consequence, most targets are located in 
Ophiuchus, Serpens, Lupus and Corona Australis. 

For any given single observation, the wavelength cover- 
age will not be complete due to the presence of saturated 
telluric absorption lines. This is particularly true for CO 
and water, for which the Earth's atmosphere absorbs in 
the same transitions as those targeted. For this reason, 
a number of observations were repeated with a cadence 
of 3 months to more than a year. The primary purpose 
was to take advantage of the Earth's velocity around the 
Sun to shift the telluric CO lines relative to those of the 
target disks. By observing targets at different epochs, 
complete line profiles could be constructed by combining 
of spectra obtained at each side of a given object's transit 
date. An example of how observations at two different 
epochs were combined to complete the spectral coverage 
of the CO lines is shown in Figure 1. This strategy also 
allows for a shallow search for variability in the astro- 

^ A correction for differential refraction was implemented in the 
observing software on November 28, 2008. 
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Figure 1 . Example of how observations at several epochs are com- 
bined to fill gaps left by saturated telluric lines. The flux spectra 
(left) and spectro-astrometry (right) shown are of VV Ser taken at 
a slit P. A. = 15°. 

metric spectra, although the shifting telluric absorption 
often make sensitive comparisons difficult. Variability at 
such time scales may be expected since the Keplerian 
time scale at the radii traced, ~ 0.1 — 1 AU is similar 
to or shorter than the cadence of observations. We can 
note that we did not detect any obvious variability in the 
sources observed during several epochs, but that we do 
not rule variability below the ~20% level. 

Significant artifacts in the astrometric spectra due to 
flat field and point spread function (PSF) effects may 
remain. The latter occurs if the PSF is not exactly rota- 
tion symmetric - and PSFs never are. Indeed, Brannigan 
et al. (2006) found such artifacts to be a common fea- 
ture of spectro-astrometric observations. The correction 
for PSF artifacts is therefore an essential calibration of 
spectro-astrometry, without which meaningful analysis is 
not possible. To effectively correct for PSF artifacts, all 
astrometric spectra were obtained with the slit oriented 
at the desired position angle (P. A.), as well as at an anti- 
symmetric P.A. + 180°. In order to ensure that the instru- 
ment was kept as stable as possible, a special CRIRES 
observing template was developed in which the grating 
and prism angle piezos were kept unchanged while the 
derotator changed the slit position angle. Furthermore, 
identical jitter patterns were used for the parallel and 
anti-parallel slit positions to maximize the reproducibil- 
ity of artifacts. The difference average between these 
two antisymmetric spectra cancel out PSF artifacts while 
preserving any real signal. The efficacy of the procedure 
is demonstrated in Figure 2. 

2.2. Data reduction 

The data were reduced using our own IDL scripts. The 
procedure includes flat-fielding, correction for the non- 
linearity of the CRIRES detector response following the 
description in the CRIRES documentation issue 86.2, co- 
adding individual nod pairs and correction for spatial dis- 
tortion. The fiux spectra were extracted using optimal 
extraction (Horne 1986), and were corrected for telluric 
absorption by division with a spectrum of an early-type 
standard star. The telluric standards were corrected for 
small airmass differences using a simple Beer law by min- 
imizing the telluric noise in a region of the spectrum rel- 
atively clear of intrinsic lines. The CRIRES grating posi- 
tion is not reproducible, so small relative shifts of order a 
few pixels in the dispersion direction between the science 
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Figure 2. Example of how parallel and anti-parallel slit positions 
are used to calibrate the astrometric spectra. Shown is the P.A.= 
55° CO observation of AS 205N. It is seen that a single spectrum is 
filled with astrometric artifacts, in this case mostly due to telluric 
O3 and CO. Also apparent is a residual low frequency artifact due 
to uncertainties in the distortion correction, which is also removed 
by the self-calibration, producing a flat astrometric flnal spectrum 
with an RMS of ~0.5mas. 

and telluric spectra were applied. Finally, in some cases 
small differences in effective resolving power between tar- 
get and standard star observation were corrected by de- 
grading the resolving power by 0.1 — 0.2 kms~^ in either 
the science or telluric spectrum. Such differences may 
arise if differences in AO correction cause differences in 
the degree to which the source fills the slit. No attempt 
was made to include an absolute flux calibration and all 
spectra were normalized to the continuum flux level. 

To further improve the signal-to- noise ratio, lines from 
similar transitions were averaged, wherever possible. In 
particular, this is important for the primary survey of the 
CO rovibrational band around 4.7 /im, where typically 8 
nearly identical (in terms of energy and collisional rates) 
lines can be averaged in a single CRIRES setting. The 
uncertainty in the spectro-astrometric signals for typical 
sources is background-limited. 

2.3. Spectro-astrometric formulation 

For the spectro-astrometry, we find that the most nu- 
merically stable and reproducible centroid (X) estimator 
as a function of line velocity v is of the form: 



X{v) = K 



[pixels] (1) 



where Xi is the spatial location of a pixel and Fi is the 
flux contained in that pixel. The centroid X(v) must 
also be corrected, using a constant factor for the fact 
that not all light will be included within the aperture. 
That is, X{v) depends on the range over which i is de- 
fined in the centroid estimator. In practice, however, K 
is small, < 1.5, and is estimated using a modeled PSF 
measured on the continuum of the spectrum. Conse- 
quently, the relative accuracy on the amplitude of the 
spectro-astrometric signal is likely < 10%, which is con- 
firmed by repeated observations of the same sources (see 
also §2). 

The uncertainty on the centroid is: 



a\X) 



where the dependency of X and F on v is left out for 
clarity. 

At this point, there is one more issue to consider. In 
our formulation, the flux is a sum of a continuum term 
and a line term: 

F^=Fc,^^FL,^. (3) 

Normally, one would be interested in determining the 
centroid offsets for the line term Fl{v) only; the presence 
of the continuum term causes X{v) to underestimate the 
true spatial offset. We call this effect continuum dilution 
(see also Pontoppidan et al. 2008). Algebraic manipula- 
tion of equations 1 and 3 shows that the relation between 
the measured and true centroid offset is: 

Xtrueiv) = Xobs(^) X (1 + Fc (v) / Fl{v)) . (4) 

It is seen that the centroid naturally diverges for Fl ^ 0. 
This makes it inconvenient to display the observed spec- 
tra in terms of Xtrue('^): and we consequently display 
all spectra in terms of Xobs('^)- The reader should thus 
be aware that the true spatial extent of the line emis- 
sion is higher by a factor 1 + Fc{v)/Fl{v). In the fol- 
lowing, we deflne a scalar amplitude of a given astro- 
metric spectrum as the maximal value of the centroid 
offset: A = Max(|X('u)|). This is a convenient model- 
independent observable that provides a measure of the 
size of the line emitting region. 

3. KEPLERIAN SOURCES 

3.1. Simple geometric models 

Keplerian sources are characterized by double-peaked 
line profiles in combination with broad astrometric spec- 
tra that display an anti-symmetric pattern at all slit posi- 
tion angles. For disks viewed at inclinations close to face- 
on, the double-peak may blend into a single peak, but if 
the astrometric spectra show an antisymmetric structure 
at all position angles, they are still considered Keplerian. 
The three disks discussed in Pontoppidan et al. (2008), 
HD 135344B, SR 21 and TW Hya, are in this category. 
Prototypical disks with clean double-peaked structure in- 
clude GQ Lup, RNO 90 and VV Ser. The astrometric 
and line fiux spectra of the Keplerian disks are shown 
in Figure 3. As shown in Pontoppidan et al. (2008) this 
structure is well explained by a simple model of a radial, 
flat distribution of gas in circular orbits around a point 
mass, and the antisymmetric structure is due to the rel- 
ative spatial displacement of red- and blue-shifted gas in 
an inclined disk. If the slit is oriented along the major 
axis of the projected disk, the maximal astrometric dis- 
placement amplitude is seen. Conversely, for an exactly 
axisymmetric and flat disk, a slit oriented along the disk 
minor axis produces no astrometric signal. Given at least 
two slit position angles, the disk position angle can be de- 
termined with a high degree of confldence. The dominant 
line emitting radii can be determined, using the maximal 
astrometric offset for a slit aligned along the disk major 
axis. Finally, the disk inclination can be determined with 
confldence if a stellar mass is assumed and vice versa. 

We use the same simple geometric model as Pontoppi- 
dan et al. (2008) to flt the data and determine the basic 
geometric parameters for the Keplerian disks. The pa- 
(2) rameters varied are the inner emitting radius, the stellar 
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Figure 3. Keplerian model fits using the simple fiat-disk model superimposed on the data. The fiux spectra are shown on top, while the 
lower panels show the spectro-astrometry at different slit position angles. The red curves show the model fits. Note that the y-axes are 
given in continuum-diluted units, and the real physical extent probed by the spectro-astrometry is larger by factors of (1 + Fq/Fl). 



mass, the disk inclination, z, and position angle, P. A. 
Figure 4 shows the resulting goodness-of-fit contours for 
these four parameters. Table 3 summarizes the best-fit 
parameters. As expected, the stellar mass and the disk 
inclination are degenerate, but in such a way that even 
an uncertain assumption of the stellar mass allows an 
accurate determination of the disk inclination. The disk 
position angles are absolutely determined, while the in- 
ner radius is somewhat dependent on the choice of the 
T{R) relation (here assumed to be a power law with ex- 
ponent q = —0.5). In lieu of the inner radius determined 
from the fit, the size of the line emitting ring, in an aver- 
aged sense, can also be estimated by using the amplitude 
A, defined above, and correcting for continuum dilution 
(see §2.3). 

3.2. A size-luminosity relation for ro-vibrational CO 

In Figure 5, the astrometric amplitudes are plotted ver- 
sus the stellar luminosities. It is seen that the Keplerian 
disks show a clear correlation across nearly four decades 
in stellar luminosity as A oc L^, with a best-fit exponent 
of a = 0.48. This is the relation expected for the radius 
of a specific equilibrium temperature as a function of 
stellar luminosity (e.g., Dullemond et al. 2001; Monnier 
& Millan-Gabet 2002). A very similar size-luminosity 
relation was found for the i^-band continuum emission 
from Herbig Ae/Be stars using the Keck interferometer 
(Monnier et al. 2005). Here, we find that the molecu- 
lar gas obeys a similar relation, but on larger scales and 
extending down to sub-solar luminosities. 

The observed points are compared to the radii of dif- 
ferent optically thin grey dust equilibrium temperatures, 
Rs = l.ly^(L,/103 L0)i/2(T,/15OO K)-'^ (Monnier 
& Millan-Gabet 2002), where eg is the ratio of the dust 
absorption efficiencies at the color temperatures of the 
incident and reemiiting radiation fields. Note that the 
dust temperature also depends on additional radiation 
sources, such as that of the surrounding disk (Dullemond 
et al. 2001), which will tend to push the radius at a spe- 
cific temperature outwards. However, this prescription 
allows a direct comparison to the analysis of Monnier 
et al. (2005). 

As is seen in Figure 5, the gas lines are indeed domi- 
nated by gas at radii well beyond the interferometrically 



measured dust sublimation radii at T^IOOO — 1500 K, 
and consistently match dust at 350 K. In the limit of a 
disk truncated at a sharp, optically thick, inner disk rim, 
the CO sizes correspond to dust at somewhat higher tem- 
peratures of ~500K. It is important to realize that some 
molecular gas may still extend inwards, as indicated by 
the best-fit inner radii and the high velocities of emis- 
sion in the line wings (Salyk et al. 2007), but the lines 
are not dominated by that component; the astrometric 
sizes measure the bulk of the gas emission. 

Of particular interest, however, is that some transition 
disks fall on the size-luminosity relation defined by the 
classical disks, specifically TW Hya and HD 135344B - 
SR 21 falls somewhat above the relation. This is consis- 
tent with the findings of Salyk et al. (2007), Pontoppidan 
et al. (2008) and Salyk et al. (2009); that the inner disk 
gas of some transition disks often has not been removed 
in the same way as the population of small dust grains. 
For these disks, it now appears that there is not even a 
difference in the radii forming the CO gas emission, and 
that the lines simply follow the expected dependence on 
luminosity. This can be interpreted as the removal of 
the small dust opacity component through the process 
of planetesimal formation, as this will preserve a signifi- 
cant column of gas in the inner evolved zone of the disk. 
The caveat is that this is a very small sample of tran- 
sition disks, and the inclusion of additional disks may 
show a greater degree of complexity, including the oper- 
ation of other mechanisms of inner disk clearing, such as 
photoevaporation or clearing by an unseen, but massive, 
(stellar) companion, both of which would tend to move 
the inner edge of the gas-disk outwards. It will be inter- 
esting to see how many transition disks, in fact, fall on 
the relation. 

3.3. Spectro-astrometry of hydrogen recombination lines 
in Keplerian sources 

Spectro-astrometry suffers from a fundamental sym- 
metry ambiguity. Because the line centroid offsets are 
measured relative to the continuum, departures from cir- 
cular symmetry in the continuum brightness distribution 
will be imprinted in the astrometric signal. For instance, 
if the continuum emission is due to a sharp inner rim in 
the dust disk (as in the models of Dullemond et al. 2001) 
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Table 1 

Log of observations 



Star 


PA 


Obs. Date Int. Time 


Spectral Range 






minutes 


/im 



LkHa 330 
LkHa 330 
LkHa 330 
CW Tau 
CW Tau 
CW Tau 
DR Tau 
DR Tau 
DR Tau 
TW Hya 
TW Hya 
HD 135344B 
HD 135344B 
HD 135344B 
GQ Lup 
GQ Lup 
GQ Lup 
GQ Lup 
GQ Lup 
GQ Lup 
HD 142527 
HD 142527 
RU Lup 
RU Lup 
RU Lup 
RU Lup 
HD 144432 
HD 144432 
HD 144432 
AS 205N 
AS 205N 
AS 205N 
AS 205N 
AS 205N 
AS 205N 
AS 205N 
AS 205N 
DoAr24E S 
DoAr24E S 
DoAr24E S 
SR 21 
SR 21 
SR 21 
RNO 90 
RNO 90 
RNO 90 
VV Ser 
VV Ser 
VV Ser 
VV Ser 
VV Ser 
S CrA N 
S CrA N 
S CrA N 
R CrA 
T CrA 
T CrA 



0° 

60° 

120° 

0° 

60° 

120° 

0° 

60° 

120° 

63° 

153° 

0° 

60° 

120° 

0° 

60° 

120° 

0° 

60° 

120° 

60° 

150° 

0° 

0° 

60° 

120° 

6° 

66° 

126° 

55° 

115° 

55° 

115° 

175° 

55° 

115° 

175° 

30° 

90° 

150° 

10° 

70° 

130° 

0° 

60° 

120° 

15° 

75° 

15° 

75° 

135° 

30° 

90° 

150° 

0° 

0° 

90° 



29/12/2008 24 
29/12/2008 24 
29/12/2008 24 
1/1/2009 24 
1/1/2009 24 
1/1/2009 24 
14/10/2007 24 
14/10/2007 24 
14/10/2007 24 
26/4/2007 40 
26/4/2007 
22/4/2007 
4/9/2007 
5/9/2007 
2/5/2008 
2/5/2008 
2/5/2008 
4/8/2008 
4/8/2008 
4/8/2008 
7/8/2008 
7/8/2008 
26/4/2007 
27/4/2008 
27/4/2008 



40 
20 
20 
20 
24 
24 
24 
24 
24 
24 
12 
12 
20 
24 
24 



27/4/2008 24 
2/8/2008 
2/8/2008 
2/8/2008 
29/8/2007 
29/8/2007 
29/8/2007 
29/8/2007 
1/9/2007 
2/5/2008 
2/5/2008 
2/5/2008 
3/9/2007 
2/9/2007 
2/9/2007 
30/8/2007 
30/8/2007 
31/8/2007 
25/4/2007 
26/4/2007 
26/4/2007 
5/9/2007 
5/9/2007 
1/5/2008 
1/5/2008 
1/5/2008 
4/9/2007 
4/9/2007 
3/9/2007 
1/9/2007 
26/4/2007 
26/4/2007 



16 
16 
16 
16 
16 
16 
16 
24 
16 
16 
16 
20 
20 
20 
32 
32 
32 
16 
16 
16 
20 
20 
32 
32 
32 
20 
20 
20 
12 
20 
20 



4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.805- 
4.805- 
4.805- 
4.660- 
4.660- 
4.645- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.639- 
4.639- 
4.660- 
4.660- 
4.660- 
4.660- 
4.639- 
4.639- 
4.639- 
4.660- 
4.660- 
2.905- 
2.905- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 
4.660- 



4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.901 
4.901 
4.901 
4.770 
4.770 
4.755 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.749 
4.749 
4.770 
4.770 
4.770 
4.770 
4.749 
4.749 
4.749 
4.770 
4.770 
2.977 
2.977 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 
4.770 




1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 
GQLUP CO v=l-0 




0.6 0.8 1.0 1.2 



HD142527 CO v=l-0 




RNO 90 CO v=l-0 




VV Ser CO v=l-0 





GQLUP CO v=l-0 




0.10 0.15 0.20 
R [AU] 



HD142527 CO v=l-0 




0.2 0.4 0.6 0.8 

R. [AU] 



HD144432 CO v=l-0 




RNO 90 CO v=l-0 




0.05 0.10 0.15 0.20 



VV Ser CO v=l-0 




2.0 2.5 3.0 3.5 4.0 4.5 5.0 

R. [AU] 



Figure 4. Goodness-of-fit surfaces, using the statistic for the 
simple Keplerian model fits to the Keplerian sources (see also Pon- 
toppidan et al. 2008). 



and the disk is viewed at some inclination, an asymme- 
try should be seen in the astrometric line spectrum when 
the slit is oriented along the disk minor axis. However, 
it is fundamentally not possible to disentangle this effect 
from an azimuthal asymmetry in the line intensity. In 
other words, CO spectro-astrometry cannot distinguish 
between structure in the spatial distribution of line inten- 
sity and the spatial distribution of continuum intensity. 
This could be remedied if there were an independent line 
tracer of the stellar location. In this section, we sug- 



gest that hydrogen recombination lines may, for many 
young stars, constitute such a tracer that may allow for 
spectro-astrometry of the continuum. 

There is a growing consensus that a dominant con- 
tributing process to H I emission from T Tauri and Her- 
big Ae stars is magnetospheric accretion (Calvet & Hart- 
mann 1992; Muzerolle et al. 1998; Bary et al. 2008), a 
shift from an original interpretation in the framework of 
stellar winds (Hartmann et al. 1990; Calvet et al. 1992; 
Grinin & Mitskevich 1991), based on observed P Cygni 
profiles of the Balmer lines in some sources. While winds 



7 



Table 2 

Disk and stellar properties 



Source 


Class^ 


line profile^ 


distance^ 


^CO(LSR) 




Sp. T. 






references^ 








DC 


km s""*^ 












LkHa 330 


trans, disk 


Keplerian 


250 


9.0 


16 


G3 


2.5 


-8.80/-8.80 


2,12 


CW Tau 


CTTS 


self-abs. 


140 


7.5 


0.8 


K3 


1.2 


-8.80/-7.99 


15,19 


DR Tau 


CTTS 


single-peak 


140 


10.7 


0.9 


K5 


1.0 


-7.5/-5.1 


16,17,18,19 


TW Hya 


trans, disk 


Keplerian 


51 


3.0 


0.23 


K7 


0.7 


-8.80/-8.80 


14 


HD 135344B 


trans, disk 


Keplerian 


84 


7.5 


8 


F3 


1.6 


-8.30/-8.30 


6 


GQ Lup 


CTTS 


Keplerian 


150 


1.0 


0.8 


K7 


0.8 


-8.00/-8.00 


11 


HD 142527 


HAeBe 


Keplerian 


198 


5.0 


69 


F6 


3.5 


-7.16/-7.16 


6 


RU Lup 


CTTS 


single-peak 


150 


3.5 


0.4 


K7 


0.7 


-7.70/-7.70 


5 


HD 144432 


HAeBe 


Keplerian 


145 


6.0 


10 


A9 


1.7 


7.07/-7.07 


6 


AS 205N 


CTTS 


single-peak 


125 


4.5 


7.1 


K5 


1.1 


-6.14/-6.14 


10 


DoAr 24E S 


CTTS 


self-abs. 


125 


3.5 


1.3 


K7-M0 


0.7 


-8.46/-8.46 


8 


SR 21 


trans, disk 


Keplerian 


125 


3.0 


15 


G2.5 


2.2 


< -8.84 


2 


RNO 90 


CTTS 


Keplerian 


125 


-1.5 


4.0 


G5 


1.5 




3 


VV Ser 


HAeBe 


Keplerian 


415 


7.0 


125 


B1-A3 


3.0 


-6.34/-6.34 


9,4 


S CrA N 


CTTS 


single-peak 


130 


2.4 


2.3 


K3 


1.5 




13 


R CrA 


HAeBe 


self-abs. 


130 


5: 


100 


B8-F5 


3.5: 


-7.12/-7.12 


1 


T CrA 


CTTS 


self-abs. 


130 


7: 


8 


F0-F5 


1.6: 


< -8.20 


7 


References. 


— [1] Bibo et al. (1992), [2] Brown et al. (2007) 


, [3] Chen et al. (1995), [4] Dzib et al. (2010), 



[5] Herczeg & Hillenbrand (2008), [6] Meijer et al. (2008), [7] Meyer & Wilking (2009), [8] Natta et al. (2006), [9] 
Pontoppidan et al. (2007a), [10] Prato et al. (2003), [11] Seperuelo Duarte et al. (2008), [12] Salyk et al. (2009), [13] 
Schegerer et al. (2009), [14] Thi et al. (2010), [15] White & Chez (2001), [16] Mora et al. (2001), [17] Muzerolle et al. 
(2003), [18] Gullbring et al. (2000), [19] Johns-Krull & Gafford (2002) 

^ Type of disk - can be classical T Tauri star (CTTS), Herbig AeBe star or transition disk 
^ Type of CO rovibrational line profile as discussed in §1.2. 

^ The distances are based on the current best estimates to the parent young clusters of the disks, many of which 
are determined by parallax measurements of known cluster members. (Dzib et al. 2010; Torres et al. 2009; Loinard 
et al. 2008), (Torres et al. 2009). The distance to Corona Australis is well determined using the orbit solution for 
the eclipsing binary TY CrA (Casey et al. 1998). One exception is HD 135344B, which has an uncertain distance of 
84-140 pc (Grady et al. 2009). 

^ The mass of the central star is estimated based on the luminosity and spectral type using the evolutionary tracks of 
Siess et al. (2000). 

^ Range of mass accretion rates found in the literature. 

^ References used for the stellar properties and the mass accretion rates. 

in Figure 6. No astrometric signals are detected to lim- 
its of 0.35 mas for RNO 90, while a tentative astromet- 
ric signal is seen at the 0.2-0.3 mas level at a P.A= 60° 
for HD 142527. The formal displacement errors in the 
two sources correspond to similar physical sizes since 
HD 142527 is at almost twice the distance of RNO 90 
(198 pc versus 125 pc). The conclusion is that the dust 
continuum emission is azimuthically symmetric and spa- 
tial asymmetries seen in the CO emission is not due to 
spatial structure of the continuum emission. 

These results are consistent with an accretion origin of 
the infrared H I lines. Alternatively, formation in a very 
compact stellar wind cannot be ruled out (a slight blue- 
shifted asymmetry is seen in the line profiles). Note that 
a potential stellar wind giving rise to H I lines should 
not be confused with the much more extended disk wind 
discussed below in the context of the CO lines. More 
importantly, in the context of the inner disk, the lack 
of strong astrometric signatures from the H I lines rule 
out a sharp, inclined inner edge structure of the 4 /im 
continuum emission. This is a result similar to that found 
with near-IR interferometry of the 2 /im continuum. It 
is interesting to note that Whelan et al. (2004) found 
highly extended (10-100 mas) spectro-astrometric signals 
of the Paschen /3 lines (tracing the same n=5 level as the 
Brackett a line used here) from a few T Tauri stars, 
including DG Tau, which is known to possess a strong 
jet. 

4. NON-KEPLERIAN (RADIAL) FLOWS 



Table 3 

Best fit model parameters for Keplerian sources 



Star [AU] P.A. 



LkhA 330 


4±1 


218±10° 


12±2° 


GQ Lup 


< 0.1 


357±10° 


65±10° 


HD 142527 


0.2±0.3 


299±3° 


20±2° 


HD 144432 


2.7di0.1 


95di3° 


25di3° 


RNO 90 


0.06di0.01 


177di3° 


37di4° 


VV Ser 


0.3±0.3 


17di4° 


65di5° 



^Assuming the stellar masses from Table 2 

likely do contribute to the optical H I lines, near-IR lines 
tracing higher energies and densities appear to be dom- 
inated by accretion flows, as indicated by a general lack 
of blue-shifted absorption (Folha & Emerson 2001). If 
the lines are indeed dominated by accretion flows, the 
astrometry is expected to constrain the H I emission to 
within a few stellar radii, corresponding to scales signif- 
icantly smaller than the disk co-rotation radii. 

We obtained spectro-astrometry along two perpendic- 
ular position angles of the Brackett a lines at 4.05 /im for 
two of the Keplerian sources: RNO 90 and HD 142527. 
This line traces somewhat lower energies than e.g. the 
Brackett 7 line at shorter wavelengths, 2.16 //m, but was 
chosen because the 4 /im continuum source is more likely 
to be comparable to that at 4.7 /im. In particular, it is 
more likely to trace dust emission, which is not neces- 
sarily the case at 2 /im, where gas opacity may dominate 
(Eisner et al. 2009). The resulting spectra are shown 
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Figure 5. Relation between the spectre- astrometric amplitude 
(continuum-corrected) and stellar luminosity for the four classes of 
astrometric spectra. The astrometric amplitudes have been cor- 
rected for continuum dilution and are normalized to a distance of 
125 pc. Top: The Keplerian profiles divided into the three tran- 
sition disks from Pontoppidan et al. (2008) and the CTTS's from 
this paper. The Keplerian disks exhibit a tight correlation with 
A oc LJ-^. Bottom: The same as the top plot, but for the single- 
peaked and self-absorbed sources. No correlation is seen. Further, 
the self-absorbed sources show much greater amplitudes than the 
Keplerian disks. 

Only some CO spectra of protoplanetary disks are 
as simple to interpret as the Keplerian disks discussed 
above. Many show a structure not consistent with a 
strictly defined Keplerian velocity field. The class of 
single-peaked line sources is characterized by a single 
peaked line spectrum, but with a broad base extending 
to ~ 50kms~^ and narrow astrometric spectra that are 
highly asymmetric about the continuum at a character- 
istic position angle. This spectral class was described 
by Najita et al. (2003). Bast et al. (2010) notes that 
the narrow central peaks in these sources do not show a 
splitting at least down to the CRIRES spectral resolution 
of ~ 3kms~^. One obvious explanation for the narrow 
peak, favored by Occam's razor, involves a Keplerian disk 
with line emission from large radii with corresponding 
low Keplerian velocities. The spectro-astrometry results 
rule out this scenario. First, if the lines are formed in a 
Keplerian disk, the central peak must be emitted from 
radii of i^peak ~ GM*(2 sin z/'u^j^jp^es)' corresponding to 
> 30 AU for i > 15°. Such extended emission should 
be directly spatially resolved by classical imaging with 
CRIRES (which has a diffraction-limited 4.7 /im spatial 
resolution of ~ 20 AU at 125 pc), yet no extended line 
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-0.1 
-0.2 




PA=0 degrees 



PA=90 degrees 



200 
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Figure 6. Spectro-astrometry of the H I Brackett a lines for 
RNO 90 and HD 142527. Because the H I hnes are highly over- 
resolved at the CRIRES resolution, the spectra have been rebinned 
to a sampling of 12kms~-'^ to maximize the signal-to-noise. The 
formal noise on the HD142527 astrometric spectrum is ~0.13 mas, 
or ~ 0.1 AU at a distance of 198 pc, including continuum dilution. 
The vertical dashed lines indicate the velocities of the CO rovibra- 
tional line centers. 

emission is observed. However, classical imaging still al- 
lows for the possibility of disks with i < 15°. This is 
where spectro-astrometry steps in by constraining the 
line emission to much smaller spatial scales (< lAU), 
regardless of inclination, thus ruling out formation in a 
Keplerian flow at the ~100<j level. 

AS 205N is the prototypical example of a source with 
single-peaked, broad winged line proflle. The astromet- 
ric spectra of the three clean single-peaked sources are 
shown in Figure 7. For these sources, one out of three 
position angles, separated by 60° shows an astrometric 
line signature centered on the line velocity, but entirely 
offset in one direction. The other two angles show the 
antisymmetric signature that might be expected from a 
Keplerian velocity fleld. The amplitude of the astromet- 
ric spectra are < 1 AU, and the astrometric line is sig- 
niflcantly narrower than the flux proflle, with the wings 
missing. 

In essence, the combination of narrow line proflles and 
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Figure 7. Flux spectra (top) and spectro-astrometry (lower panels) for the single-peaked line sources. The red curves show possible wind 
models (see Table 4) 



small spatial extent indicates sub-Keplerian motions. We 
explore a model that explains the qualitative properties 
of the single-peaked class of CO line spectra in §5. 

4.1. H2O and OH in AS 205N 

A number of other single-peaked line disks are known 
to show strong emission from water vapor at 3 /im as well 
as throughout the mid-infrared range. Salyk et al. (2008) 
demonstrated the existence of lines due to gaseous H2O 
and OH in the 3 /im hot band for AS 205N. The lines 
have low contrast relative to the continuum (5-10%). 
However, CRIRES is significantly more sensitive at 3 /im 
than at 4.7 /im, so a detection of an astrometric signal 
from the water lines may be possible. AS 205N was ob- 
served in August 2007 with 3 different position angles to 
match the CO observation. No astrometric signals were 
detected either for H2O or for OH. The upper Hmits are 
in both cases 0.2 x (1 + Fc/F/)mas, corresponding to 
< 0.5 AU at a distance of 125 pc, when correcting for con- 
tinuum dilution, or an emitting area of < 0.8 AU^. This 
is consistent with a measured emitting area of 0.4 AU^ 
for the 3 /im water lines found by Salyk et al. (2008), 
and is marginally smaller than the radial extent of the 
CO emission (0.7 AU). Since the excitation energies of 
the 3 /im water and OH lines are higher than those of 
CO, it is likely that the smaller extent can be explained 
by an origin in somewhat warmer gas. 

4.2. Self-absorbed sources 

Finally, self-absorbed sources have strong absorption 
components superposed on broad emission lines and are 
consequently more complex. While the gas forming the 
absorption lines in some sources may be completely un- 
related to the disk, it is discussed in section 5 how some 
self-absorbed sources may be a different representation of 
the peaky line sources, namely those viewed at a higher 
inclination angle. 



The self-absorbed sources tend to show the highest am- 
plitude astrometric spectra of the survey. The physical 
difference between the single-peaked and self-absorbed 
sources and the Keplerian disks is illustrated by their as- 
trometric amplitudes in Figure 5. These sources do not 
fall along a neat correlation with the stellar luminosity 
as do the Keplerian disks. T CrA and R CrA in partic- 
ular show astrometric amplitudes as high as 20-30 AU, 
indicating that the absorbing gas is structured on much 
larger scales than the emitting gas. These sources are 
likely to be younger and more embedded than the re- 
maining sample, which may impact the observables. 

5. A NEW UNIFIED DISK+WIND MODEL 
5.1. Parametrization 

A parametrized disk wind model is used to test whether 
an idealized wind structure can produce a plausible 
match to the phenomenology of the single-peaked line 
profiles and associated asymmetric astrometric spectra. 
It is based on models for UV resonance and hydrogen line 
emission from accretion disk winds (e.g., Knigge et al. 
1995; Kurosawa et al. 2006), and represents a computa- 
tionally convenient structure inspired by numerical re- 
sults from magneto-hydrodynamical (MHD) simulations 
of magnetized disk winds (Konigl & Pudritz 2000). 

The basic requirement of the observed spectro- 
astrometry is that the fine forming gas must be orbiting 
the central star with strongly sub-Keplerian azimuthal 
velocities in order to produce the single peak without 
requiring that the emission is extended at the spatial 
resolution of CRIRES (~0V15). A wide angle wind pro- 
vides a convenient physical way of accomplishing this 
through simple conservation of angular momentum - as 
a gas parcel is forced outwards due to the wind pressure, 
the azimuthal velocity decreases linearly with radius, in 
comparison with the underlying Keplerian disk in which 
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Figure 8. Spectro-astrometry of water and OH lines from the AS 
205N disk. 



the velocity experiences a shallower decrease as R~^^^. 
This generates gas above the disk that is supported by 
wind pressure and orbits at low azimuthal velocities. 

Following Kurosawa et al. (2006), the wind is con- 
structed as a set of linear streamlines with a locus below 
the central star at a distance d in units of . This gen- 
erates a conical wind with no flow along the disk axis. 
Briefly, the wind is accelerated along the field lines as: 

'^stream = + (/ * t'esc ~ Cg) X (1 — Agc^le/^ + ^scale)^ 7 (5) 

where / is the coordinate along the stream line, Cg is the 
sound speed, Vqsc is the asymptotic velocity at the end 
of the stream line and A^caie is the scale of the accel- 
eration region of the wind. /3 is the wind acceleration 
parameter. Requiring angular momentum conservation, 
the azimuthal velocity component is: 



Vc/y = VKepler(^o) X 



Rjlo) 
R ' 



(6) 



conservation: 

Pwind 



t{R) 



^stream I COS^| 



S cosS 



(7) 



Here, il(i^) oc R~p is the local mass- loss rate, S is the 
angle between the stream line and the disk normal and 
S is the distance to the wind locus. The exponent of the 
local mass loss rate is taken to be p = 7/2 (Krasnopol- 
sky et al. 2003). The total wind mass loss rate can be 
calculated by integrating over the disk and multiplying 
by two to include the opposite surface: 



M = 2 / 



27iRt{R)dR 



(8) 



where R is the radial disk coordinate. 

The density of the wind is calculated assuming mass 



The raytracer RADLite (Pontoppidan et al. 2009) is 
used to render model lines and spectro-astrometry for the 
wind models, based on a generic model of a flared proto- 
planetary disk, and assuming level populations in LTE. 
Specifically, the temperature structure is assumed to be 
in equilibrium with the stellar radiation field and dom- 
inated by dust heating/cooling. In reality, the heating 
of the wind is likely to be dominated by photo-electric 
heating similar to the heating of the disk atmosphere 
(Kamp & Dullemond 2004; Jonkheid et al. 2004; Gorti 
& Hollenbach 2004; Dullemond et al. 2007) or, perhaps, 
ambipolar diffusion (Safier 1993). The cooling may be 
dominated by adiabatic expansion and molecular cool- 
ing (e.g., partly via the observed CO and H2O lines). 
However, we restrict ourselves to qualitative models in 
this paper (see also §6.5), since a detailed and appro- 
priate physical treatment of the thermal wind structure 
required to match the observations will be likely be a 
significant study in its own right. 

5.2. Properties of the wind model ohservahles 

Figure 10 illustrates the wind geometry and com- 
pares the observables generated using the wind model 
for the spectro-astrometry of AS 205N. The total mass- 
loss rate is 9 x 10~^ Mq yr~^, assuming a CO abundance 
of 5 X 10-^ relative to H. This mass-loss rate is con- 
sistent with a highly accreting CTTS. Indeed, AS 205N 
is accreting at a rate of 7 x lO^^M© yr~^ (Eisner et al. 
2005). There is no treatment of the chemistry of the 
wind in the present model, so if the abundance of CO is 
lower than that assumed, the mass-loss rate will be cor- 
respondingly higher. However, the numbers appear to 
be physically reasonable; the ratio of mass-loss to mass- 
accretion rates is Mwind/^acc ^ 1.5%, and as this ra- 
tio is expected to be as high as 10% (Konigl & Pudritz 
2000), the CO abundance in the wind could in fact be 
lower than assumed. There are many free parameters in 
the wind model, so that shown represents one possibility 
that reproduces the main characteristics of the observed 
astrometric spectra. A full parameter study, beyond the 
scope of this paper, may reveal significant degeneracies. 
That said, some parameters are well constrained, includ- 
ing the position angle of the system. General properties 
are also locus distances d < 10 i?* and inclinations of 
i < 30°. The velocity field used for AS 205N is shown in 
Figure 11. 

It is seen that the wind model can explain the key 
properties of the combined line profile and astrometric 



11 



S -2 



0° 



J- 



60° 



120° 






G "-4 
o 0.2 



G 
o 

i -2 







■ 


1 




0° 









-40 -20 20 40 
V^, [km s-'] 



-40 -20 20 40 
V [km s"'] 



-60-40-20 20 40 60 
V^, [km s-'] 



-40 -20 20 40 
V^, [km s-'] 



Figure 9. Spectro-astrometry for the self-absorbed sources. Note that while a model of DoAr24E S exists that fits the spectro-astrometry, 
the fiux line profile is not well matched. One explanation for this is that there is an additional compact component to the line. No attempt 
is made to fit models to T CrA, 



spectra: namely it readily produces single-peaked lines 
with broad wings. The wings are dominated by the in- 
nermost, Keplerian disk, while the single peak is gener- 
ated by the sub-Keplerian gas at a few AU. The presence 
of a strongly sub-Keplerian component requires a wide- 
angle, non-collimated wind, in the model parametrized 
as a small star-to-locus distance, d. Increasing d results 
in a double-peaked line profile that would be resolvable 
with CRIRES. The astrometric spectra become highly 
asymmetric at P.A.'s close to the disk minor axis. It is 
important to note that this wind geometry does not nec- 
essarily produce strongly blue-shifted lines - which might 
have been expected. The uncollimated winds sees much 
of the line emission coming from gas moving on trajec- 
tories nearly parallel to the disk surface. For sources 
viewed at low inchnation angles, this generates low ra- 
dial velocities and prevents large blue-shifts, as required 
by the data. The addition of a significant amount of gas 
at high altitudes where the optical depths toward the 
central star are low, results in higher column densities of 
warm CO, leading to brighter lines with higher line-to- 
continuum ratios, a property of this class of CO spectra 
that has been observed (Bast et al. 2010). 

Table 4 summarizes the RADLite wind model parame- 
ters for the wind-dominated disks. Since we do not fully 
explore the parameter space, the entries in the table rep- 
resent a possible model, but with a caution that signif- 
icant degeneracies may exist. Further, the use of LTE 
level populations and the assumption of coupled dust 
and gas temperatures almost certainly introduce inac- 
curacies in the implied wind structures. Indeed, it was 
typically necessary to increase the luminosity of the cen- 
tral source to values well above the known stellar prop- 
erties to match the amplitude of the astrometric signal, 
especially in the case of RU Lup, the reason being that 
the gas temperatures are significantly higher than what 



can be explained by pure coupling to dust. 

5.3. The wind model at higher inclinations as a model 
for the self-absorbed disks 

As can be seen in Figure 11, there is no wind com- 
ponent along the disk axis due to the finite size of the 
inner disk rim. Because of this, a disk wind viewed close 
to face-on, as is likely the case for AS 205N, allows a 
free view of the continuum emission from the innermost 
disk, consistent with the low extinction toward the cen- 
tral star in many of our targets. However, if the incli- 
nation were higher, say ~ 45°, the entire wind column 
will be located between the disk continuum emission and 
the observer. This leads to the formation of a strong 
line absorption component, as illustrated in Figure 12. 
The Figure shows the line profiles of the prototypical AS 
205N model, but viewed at higher inclination angles. A 
deep absorption line forms at inclinations higher than 
^ 40° for this particular model, a result that can be 
compared directly to sources exhibiting a self-absorbed 
line, such as CW Tau. The line profile of this source is 
qualitatively reproduced by the AS 205N model viewed 
at an inclination angle of 55°. The depth of the absorp- 
tion line is smaller in the model, but can easily be deep- 
ened by increasing the mass-loss rate. We note that the 
small spectro-astrometric sample of self-absorbed source 
includes sources such as R CrA and T CrA are clearly 
younger and more embedded than the remaining disk 
sample, and their deeper absorption components are con- 
sistent with mass-loss rates 2-3 orders of magnitudes 
higher than that of the AS 205N model. However, there 
are more bona-fide isolated disks in our CRIRES sur- 
vey showing strong, warm CO absorption that were not 
included as spectro-astrometric targets, but with prop- 
erties otherwise resembling CW Tau. 

The four self-absorbed disks are not modeled in detail 
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Figure 10. Sketch of a disk wind model that quahtatively reproduces the broad single-peaked CO line spectra and CRIRES spectro- 
astrometry from highly accreting T Tauri stars - in this example, AS 205N, a T Tauri star in Ophiuchus. The left panels show the observed 
CO rovibrational (v = 1 — around 4.7 /im) line spectrum and spectro-astrometry at 3 different slit position angles. The middle panels 
show two 2D radiative transfer models: one for a purely Keplerian disk, and one for a Keplerian disk in which a wide-angle molecular wind 
is launched from the disk surface. The right panel contains a sketch of the basic wind geometry. 



Table 4 

Model parameters for wind-dominated sources 



Star 


P.A.^ 


ib 


^model/^**^ 


dd 


Mq yr-1 


DR Tau 


0°: 


-9° 


1.0 


1 


410-^ 


RU Lup 


80° 


35° 


90 


8 


410-9 


AS 205N 


235° 


20° 


1.8 


5 


910-9 


DoAr24E S 


235° 


20° 


1.8 


5 


910-9 


S CrA 


15° 


10° 


1.3 


4 


110-8 



^Position angle of the disk major axis. 

^Inclination of the disk rotation axis. A negative value of the 
inclination corresponds to the north pole of the disk facing the 
observer. 

^ Ratio between the luminosity of the central source used in the 
model and the actual stellar luminosity. A high ratio indicates that 
the gas is heated to temperatures significantly higher than those 
of the dust. 

^Distance between the locus of the wind and the central star (in 
units of stellar radius). d=0 corresponds to a spherical wind. 
^Mass loss rate of the wind assuming a::(CO) = 5 10-^. 

for several reasons. We found it difficult to reproduce 
the very high amplitudes of the astrometric offsets for 
the absorption component for CW Tau and T CrA, likely 
indicating the presence of extended continuum emission 
structure. DoAr 24E S has several absorption compo- 
nents that also cannot be modeled exactly with the wind 
model without the additional parameters. However, we 
do still present a best effort model of DoAr 24E S in Fig- 
ure 9 to illustrate the required complexity of this star. 



6. DISCUSSION 

The origin of CO rovibrational emission lines from pro- 
toplanetary disks has long been a puzzle. While it has 
been clear that the lines are formed close to the star, 
given the high temperatures and densities required to 
excite the lines, the great variation of the line shapes 
and excitation temperatures - rotational and vibrational 
- have evaded a unified explanation. Specifically, the 
lines cannot be explained solely by thermal emission from 
simple Keplerian disks, except in a few cases. This was 
already noted by Najita et al. (2003), but an unambigu- 
ous identification of the additional (to a Keplerian disk) 
radial flow pattern could not be determined. Now, with 
the addition of spatial information, as offered by high res- 
olution spectro-astrometry, we can present a more com- 
prehensive picture of the dynamics of molecular gas on 
AU-scales in protoplanetary disks. 

6.1. Size of the line emitting regions 

In §3 we show that sources with double-peaked (Ke- 
plerian) profiles obey a size-luminosity relation with a 
power-law index of a = 0.5, as expected for emission 
truncated at a radius of constant temperature. Criti- 
cal ingredients in identifying this relation are the abil- 
ity to distinguish Keplerian disk-dominated lines from 
the wind-dominated and self-absorbed lines, as well as 
the use of the direct size measure offered by spectro- 
astrometry. Specifically, the astrometric amplitude is not 
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Figure 11. Radial and azimuthal velocities of the gas in the 
parametrized disk+wind model. For this set of parameters, which 
reproduce the observables in, e.g., AS 205N, the wind is seen to be 
strongly sub-Keplerian beyond 0.5 AU. Note that the finite radius 
of the inner rim of the disk and a locus distance > prevents the 
wind from filling in all polar angles. In this case, the wind only 
exists at polar angles of ^ > 50° (as measured from the disk axis). 
For angles closer to the disk axis, the velocity field is still Keple- 
rian, although the gas density here is far too low to generate any 
line emission. 

dependent on disk inclination, which otherwise makes it 
difficult to use the line width as an accurate proxy for 
size. 

Another conclusion that can be drawn from the mea- 
sured sizes of the emitting region from Keplerian sources 
is that a line origin in the accretion flow can be ruled 
out. Accretion flows are expected to generated double- 
peaked CO line profiles (Najita et al. 2003). Funnel 
(along magnetic field lines) accretion flows are expected 
to be launched from radii near, or within, the corota- 
tion radius (Shu et al. 1994). This is, for T Tauri stars, 
located inside the dust sublimation radius at 2-10 stel- 
lar radii, and thus well within the observed CO radii 
of ~100i?*. The Bra lines as observed in RNO 90 and 
HD 142527 have astrometric displacements much smaller 
than those of CO and evidently traces a different gas 
component, which may still be an accretion flow. Note 
that because the H I lines trace a much wider range of 
infall velocities than the CO -HI line widths are 100- 
200kms~^ - they are less affected by Keplerian rotation 
at the bottom of the flows than the narrower CO lines 
(Muzerolle et al. 1998). Consequently, they may not dis- 
play a double-peaked profile even if formed as part of an 
accretion flow. 

Is it expected that CO follows a relation that was devel- 
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Figure 12. Effect of inclination on the line profile in the disk 
wind model. 

oped for dust subhmation? One possibility is that dust 
shielding plays a significant role in the survival of CO in 
the inner regions of disks. In this case, CO would only 
be able to survive at radii at which the radial column 
of dust is sufficiently high. This is consistent with CO 
existing at radii larger than those of the inner dust rim. 
However, CO is also know to efficiently self-shield, even 
at low column densities of N{H2) ~ lO^^cm"^ (Visser 
et al. 2009), although uncertainties are large at temper- 
atures higher than a few 100 K. 

6.2. Non- Keplerian motions and a unifying model for 
CO rovibrational lines 

We propose that most, if not all, disks have a molecular 
wind velocity component in addition to pure Keplerian 
rotation - one that is effectively traced by ro- vibrational 
CO lines, presumably in addition to lines from a wide 
range of other molecules with strong infrared emission 
bands. Najita et al. (2003) argues against a wind based 
on the fact that the lines are not blue-shifted and that 
very young (stage I) high accretion rate sources some- 
times do not show rovibrational CO emission lines. How- 
ever, our wind model shows that a strong blueshift of the 
line emission is not necessary for the uncollimated slow 
disk winds constrained by angular momentum conser- 
vation. While we do not specifically address the lack of 
strong CO lines toward some embedded young stellar ob- 
jects, strong accretion may heat the disk mid-planes to a 
point where the temperature inversion in the disk surface 
required for line emission is no longer possible. Further- 
more, we do find that the contribution of outfiowing gas 
to the line emission is highly variable, possibly in rela- 
tion to the level of accretion activity onto the central 
star. While we have not sought to determine the wind 
launching mechanism, the data do suggest that the wind 
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is uncollimated and slow (I'ywindl < '^Kepler)- It is impor- 
tant to stress that the wind should not be seen as entirely 
separate from the Keplerian disk, but rather as a mod- 
ification. For instance, the gas does not suddenly stop 
its Keplerian rotation even as it attains a radial velocity 
component. For low wind accelerations, it may be diffi- 
cult to determine that there is a difference at all from a 
pure Keplerian velocity field. As the wind is launched, 
the gas maintains its angular momentum and slows down 
so much as it is pushed outwards that it may never escape 
the system, but re- accrete onto the disk at larger radii. 
The mass-loss rates implied by our wind models are large 
enough that the inner disk gas can be completely redis- 
tributed during the lifetime of the disk. This may limit 
the time available for planet formation, in line with ideas 
generally associated with photo-evaporative disk winds 
(Alexander et al. 2006). 

6.3. Wind launching mechanism 

We will not discuss in detail how the molecular wind 
is launched. However, the requirement that it is uncol- 
Hmated and slow likely places significant physical con- 
straints on the launching mechanism. For instance, a 
magnetic centrifugal wind might be expected to generate 
winds that are much too fast as the gas in this case will 
be locked to the stellar rotation, known to be > t'Kepier 
beyond the inner rim of the disk. A thermal wind is much 
slower, but may require high ionization fractions and gas 
temperatures (Tgas ~ 10, 000 K), which could be inconsis- 
tent with the presence of abundant molecules with rota- 
tional temperatures of ~ 1, 000 K. The photo-evaporative 
wind models by Alexander et al. (2006) predict wind 
launch speeds that may even be slightly subsonic (5- 
lOkms""^), consistent with our observational finding. It 
can be noted that the evaporative winds predict a some- 
what shallower exponent p for the drop-off of the mass- 
loss rate i'{w) with radius than the p = 7/2 used here. 
Experimenting with our model indicates that good fits 
may also be found with p = 5/2, although not as readily. 
Also noteworthy is that the models by Alexander et al. 
(2006) predict significant blueshifts of wind-dominated 
[Ne II] line peaks (5-10 km s"^ Alexander 2008), which 
does not appear in the CO data, at least compared to the 
ambient cloud velocity traced by cold absorption compo- 
nents. 

6.4. Implications for radial mixing 

The observational implication of slow gas flows above 
the canonical "warm molecular layer" of the disk may 
have important implications for the transport of material 
in the disk. Specifically, the gas may never be acceler- 
ated to velocities allowing it to escape the disk, allowing 
it to fall back onto the disk at larger radii. For flow rates 
of lO~^M0yr~^, a wind lasting a few Myr will clearly 
be able to re-distribute a significant fraction, if not all, 
of the inner disk gas not accreted onto the central star. 
Similar ideas where driving the development of early disk 
wind models, such as the X-wind (Shu et al. 1994) to ex- 
plain the redistribution of solids required to produce the 
crystalline dust grains in cometary material. It may be 
important to note that if the outflowing gas falls back 
onto the disk, there will be a significant azimuthal ve- 
locity differential as sub-Keplerian material impacts the 



Keplerian disk below. Judging from Figure 11, the differ- 
ential may be as high as 10-15 km s~^ if the fall happens 
at 10 AU, but less if the fall happens farther out in the 
disk. This is sufficient to generate shocks that may be 
observable. 

6.5. Caveats, unknowns and new questions 

Clearly, the parametrized wind model is not based on 
a detailed treatment of the underlying hydrodynamics 
and radiative transfer. However, it does represent a phe- 
nomenological structure required to match clearly de- 
fined observables. It is therefore important to consider 
the circumstances under which the physical and chemical 
structure of the wind may differ from the simple model. 

One property of the wind model that does not match 
the observations is the rotational temperatures of the 
lines. The model predicts lower temperatures than ob- 
served. However, this can likely be explained by the 
use of LTE level populations in thermal balance with 
the dust. The gas that forms the wind will be exposed 
to a strong IR and UV radiation field from the central 
star and innermost disk, and is therefore subject to flu- 
orescence pumping as well as collisional excitation by a 
gas that is heated and likely partly ionized by photo- 
electrons. Balancing this are the strong cooling terms 
offered by the rovibrational molecular emission, in partic- 
ular that of CO and water. While previous wind models 
(e.g., Safier 1993) generally predict very high wind tem- 
peratures (10,000 K) and fast winds (100 km s""*^), the 
existence of significant molecular coolants may, in part 
of the wind flow, maintain the required low tempera- 
tures and velocities. In particular the temperatures may 
be closer to those given by the assumption of equilib- 
rium with the stellar radiation field, as assumed in this 
paper, than the high temperatures implied by a purely 
atomic gas. A future detailed heating-cooling balance 
calculation for the wind model is clearly an important 
next step. 

More detailed hydro-dynamic modeling is essential to 
fully assess the implications of slow disk winds. What is 
the wind launching mechanism? Does the wind eventu- 
ally fall back onto the disk and at which radii? How can 
the wind remain, at least in part, in a molecular form 
as the material is lofted to altitudes where dust shielding 
low and the gas is exposed to a harsh UV radiation field? 
While a model of the chemistry of the molecular wind is 
a study in its own right, it can be noted that the CO col- 
umn densities required to generate the deep absorption 
Hues seen in the self-absorbed (more inclined) sources are 
likely high enough to self-shield. Likewise, the presence 
of water in the flow, at least in AS 205N - see Figure 8 - 
will provide additional shielding against UV photons for 
a range of molecular species (Bethell & Bergin 2009). 

7. CONCLUSION 

Using spectro-astrometry to image molecular gas at 
0.1-10 AU in a sample of protoplanetary disks has re- 
vealed an intriguing range of structure. Some disks ap- 
pear to be dominated by gas in Keplerian rotation about 
the central star, as expected, while others have a signif- 
icant slow radial velocity component consistent with a 
wide-angle disk wind. The basic observational evidence 
for a slow wind is the clear presence of low velocity gas 
5kms~^) within a few AU from the central star. It 
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is difficult to generate low velocity gas deep in a poten- 
tial well, but one simple way of doing so is via angular 
momentum conservation of an expanding, but initially 
Keplerian, flow. While it has long been known that fast 
atomic winds were common in T Tauri stars, we now 
find that there is a significant molecular component as 
well. Further, it appears that the observed molecular 
wind must be launched from the disk surface. It is not 
clear whether the gas in the wind reaches escape veloc- 
ities, and may therefore be re- accreted onto the disk at 
larger radii. The wind reproducing the CO line spectro- 
astrometry is much slower than that predicted by X- 
wind theory, which reaches terminal velocities of several 
lOOkms"^ (Shu et al. 1994), but is a much better match 
to photo-evaporative flows that have poloidal velocities 
similar to the sound speed (e.g., Alexander et al. 2006). 
The existence of disk winds with high mass-loss rates 
have significant implications for the availability of mate- 
rial for planet formation in the PFZs of protoplanetary 
disks, and may limit the time available for planet for- 
mation. The potential large scale cycling of inner disk 
material as suggested by the low velocity of the winds will 
also influence the chemistry of planet-forming material, 
for instance by exposing a large fraction of the disk mass 
to the strong UV fields at high disk elevations. Future 
work will include the development of a model for how 
slow molecular winds are launched, as well as a chemical 
model of the wind that can explain the survivability of, 
at least, CO, H2O and OH at low optical depths and at 
high elevations above the disk. 
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APPENDIX 
NOTES ON INDIVIDUAL SOURCES 

Keplerian disks 

LkHa 330 — This transitional disks is one of the "cold disks" imaged at submillimeter wavelengths by Brown et al. 
(2009). The spectro-astrometry reveals a significantly smaller inclination angle (12°) than that suggested by the 
ellipticity of the 850 /im continuum image (42°). The submillimeter image also shows that dust has been depleted 
within ~50AU, while we find that molecular gas persists at 4 AU, but not closer to the star. In many ways, this disk 
appears to be similar to that of another transition disk, SR 21, described in Pontoppidan et al. (2008). 

GQ Lup — This star is probably best known for having a substellar companion, GQ Lup b (Neuhauser et al. 2005), 
currently located due west of GQ Lup. The spectro-astrometry of the GQ Lup disk shows that it is oriented along the 
N-S axis and has a relatively high inclination of 65 =b 10°. It is not known whether GQ Lup b orbits in the disk plane, 
but if it does, the physical separation will be higher than the projected separation of 0'!7 by a factor 2.4, bringing the 
physical separation to 240 AU. Further, the prediction is that any orbital motion of the companion will be along the 
N-S axis. It can be noted that the stellar inclination of GQ Lup has been found to be 27 ± 5° (Broeg et al. 2007) - 
very different than that of the disk, raising the question of whether this difference is due to interactions between the 
companion and the disk. 
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HD 142527 — The disk around this borderhne Herbig Ae/Be - T Tauri star (Spectral Type = F6) has been imaged in 
the continuum at near- to mid-infrared wavelengths. At mid-infrared wavelengths (18-25 /im), the outer disk exhibits 
a strongly asymmetric structure, with the eastern side being much brighter than the western (Fujiwara et al. 2006). 
These authors suggested that the bright eastern arc being an illuminated rim of a disk inclined along a north-south 
axis. The spectro-astrometry, however, demonstrate that the actual axis of the (inner) disk is at a P.A.=293°. Near- 
infrared, scattered-light imaging by Fukagawa et al. (2006) find that the disk is elliptical along a major axis in the 
NW-SE direction, more consistent with the spectro-astrometric P.A. Further, at NIR wavelengths, the south-western 
part of the disk is brighter, indicating that this is the closer part of the disk, exhibiting strong forward scattering. The 
inclination angle of 25° matches both the mid-infrared and near-infrared imaging. 

HD 144432— The polarization axis of the disk around HD 144432 was found to be 114°(Rodrigues et al. 2009). 
Monnier et al. (2006) determined the size of the H-band source (tracing an inner rim of the disk) to 0.5 AU using 
closure-phase interferometry, or slightly smaller than the CO emitting region. 

RNO 90 — While this G5 T Tauri star has featured prominently in studies of mid-infrared spectroscopic disk properties 
(Kessler-Silacci et al. 2006; Pontoppidan et al. 2010), little is known about the geometry of the disk. Using CO spectro- 
astrometry, it can be determined that the inclination is 37 zb 5° and that the major axis of the disk is oriented along 
the N-S axis with a P.A. = 177 ±3°. There is no apparent inner gap in the CO emission, which can be traced to 
within 0.1 AU. 

VV Ser — This Herbig Ae/Be star is known to have a prominent disk shadow (Hodapp et al. 2004; Pontoppidan & 
Dullemond 2005). The disk orientation, P.A. and inclination was determined by Pontoppidan et al. (2007a, b) using 
the disk shadow in combination with near-infrared interferometric visibilities. The P.A. of Pontoppidan et al. (2007b) 
of 15 ± 5 can be compared to the spectro-astrometric P.A. of 17 ± 5°, while the disk shadow inclination of 70 ± 5° can 
be compared to the astrometric inclination of 65 ± 5°. It is important to note that the likely distance to the Serpens 
star forming cloud, of which VV Ser is a member, has been significantly revised from 260 pc (Straizys et al. 1996) to 
410 pc (Dzib et al. 2010). The radius of the K-band continuum emission determined by Pontoppidan et al. (2007a) is 
then revised to 1.1 AU, in comparison to a CO emitting radius of 3.4 zb 0.4 AU. 

Single-peaked line sources 

DR Tau — DR Tau is a well-known high activity T Tauri star with a high degree of optical veiling (Hessman & Guenther 
1997) and spectral variability (Smith et al. 1997; Alencar et al. 2001). Near- infrared interferometry has provided an 
uncertain measure of the disk P.A.=160±55° (Akeson et al. 2005), consistent with our spectro-astrometry. Andrews & 
Wilhams (2007) finds a disk P.A.= 170 ± 8°, while Isella et al. (2009) finds a P.A.=98 and an inclination of 37°, using 
high resolution (sub) millimeter interferometry. Based on stellar vsini measurements coupled with periodic variability, 
MuzeroUe et al. (2003) finds a high stellar inclination of 69°. In summary, the basic geometry of the DR Tau disk is 
still uncertain, but the spectro-astrometry indicates an uncertain P.A. close to 0°and a low inclination. 

RU Lup — RU Lup is a well-known classical T Tauri star with strong accretion signatures. Gahm et al. (2005) 
suggested that the star may have a close stellar companion due to the presence of periodical radial velocity changes, 
although the presence of starspots leading to the observed radial velocity signal is not ruled out (see discussion in 
Herczeg et al. (2005)). If the variability is due to starspots, a measurement of vsini = 9 ± lkms~^ (Stempels & 
Piskunov 2002) implies a nearly face-on orientation of the star, and presumably the disk as well, of i = 24° (Herczeg 
et al. 2005). Takami et al. (2001) presented spectro-astrometry of the Ha line and found evidence for a jet with a 
P. A.- 45°. The CO P.A. of 80° is significantly different. 

AS 205N — This star, also known as AS 205A, the primary of the close binary (separation 1'.'3) system known for 
exhibiting bright molecular emission from water and organics throughout the mid-infrared (Salyk et al. 2008). Since 
the CRIRES PSF has a width of ~0'.'12, there is no contamination from the secondary. The central star has a late 
spectral type (K5) relative to its luminosity (41/©), resulting in a very young apparent age of 5 10^ yr (Andrews et al. 
2009). A comparison to the disk geometry of Andrews et al. (2009), derived with high spatial resolution submillimeter 
interferometric continuum and line imaging with the SubMillimeter Array (SMA), is instructive. They find, using the 
dust continuum contours, an inclination of ~ 25°, which is consistent with the wind model requiring the disk to be 
nearly face-on. Their P.A. is different from that derived from the spectro-astrometry. While continuum isocontour 
fits can be uncertain for nearly face-on disks, spatially resolved SMA CO (3-2) images clearly show rotation consistent 
with a disk at P.A.~165°. In contrast, the spectro-astrometry requires a P.A.~ 235° of the AS 205N disk major 
axis. We suggest that the discrepancy could be explained if the rotation pattern seen in the rotational CO line is 
due to contamination from AS 205S; a scenario which is explicitly not ruled out by Andrews et al. (2009). If so, 
this implies an interesting configuration in which the orbit of the AS 205 binary is not coplanar with the AS 205N 
disk. Non- coplanar ity has been observed as a general property of binaries with disks, although there is probably some 
alignment for disks around stars with separations between 200-1000 AU (Jensen et al. 2004). 

S CrA N — S CrA is a 1.3'' binary. Our spectro-astrometry is of the northern component, the primary at near-infrared 
wavelengths, but the secondary in the visible. There is, to our knowledge, no other measurements of the disk geometry 
of this source in the literature. 
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Self-absorbed sources 

CW Tau — CW Tau has an optical jet with a P.A.= 155° and indination of 41° (Gomez de Castro 1993; Hartigan 
et al. 2004; Coffey et al. 2008), quahtatively consistent with the CO spectro-astrometry. 

DoAr 24E S — This is a 2'' binary, with the south-eastern component being an InfraRed Companion (IRC) (Chelh 
et al. 1988). Our spectro-astrometry is of the secondary (at K-band) southern component - this is the brighter, by 
about 0.7 magnitudes, component in the M-band. The secondary disk has P.A. of ll°as determined using K-band 
polarization (Jensen et al. 2004). The secondary is itself a possible close, equal brightness, binary with separation 
5.6 mas and P.A.=44°(Koresko 2002), based on speckle imaging. The presence of a close binary may explain the 
need for two wind components components in our spectro-astrometry, separated by about 15kms~^ at the time of 
observation. It should be noted that a match with a binary is not unambiguous, and other models of extended emission 
along a P.A. =44° are also possible. 

T CrA — Bailey (1998); Takami et al. (2003) reported a binary companion, using spectro-astrometry of Ha, with a 
separation of > 0V14 and P.A. =278°. We do not detect a companion at 4.7 /im at this P.A. and separation with an 
upper hmit of the flux ratio of 6. Speckle imaging at K-band also does not detect the companion (Ghez et al. 1997; 
Kohler et al. 2008). It is possible that the companion is significantly bluer than the primary, and we do not expect it 
to influence our spectro-astrometry. 



